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Synchronization is shown to occur in spatially extended systems under the effect 
of additive spatio-temporal noise. In analogy to low dimensional systems, synchro- 
nized states are observable only if the maximum Lyapunov exponent A is negative. 
However, a sufficiently high noise level can lead, in map with finite domain of 
definition, to nonlinear propagation of information, even in non chaotic systems. 
In this latter case the transition to synchronization is ruled by a new ingredient : 
the propagation velocity of information Vp. As a general statement, we can affirm 
that if Vf is finite the time needed to achieve a synchronized trajectory grows ex- 
ponentially with the system size L, while it increases logarithmically with L when, 
for sufficiently large noise amplitude, Vf = . 



1 Introduction 

Noise-driven phenomena like diffusion and growth processes, front propaga- 
tion and interface dynamics are quite common and extensively studied prob- 
lems in equilibrium and non-equilibrium statistical mechanics El. Usually, they 
are described in terms of stochastic PDE's, or probabilistic automata, -azhere 
some random process is introduced as a model of white thermal noise an. 

More recently, noise has been shown to be able also to control chaotic 
fluctuations in low-dimensional dynamical systems, yielding synchronization 
of replicae driven lw n the same stochastic process. This problem has been 
widely investigated IJulIlJ considering two chaotic maps, e.g. the logistic one 
or the Lorenz system, coupled through the same additive noise. For sufficiently 
large noise amplitude, the invariant measure of both maps can be modified 
in such a way to favour the stabilization on the same stochastic orbit. Far 
from misterious, this effect can be explained and quantitatively predicted by 
observing that, for some noise amplitude, the Lyapunov exponent associated 
to the dynamics becomes negative, as brillantly argued by Pilovsky 0. Clearly, 
this can happen only if the map dynamics has expanding and contracting 
regions, so that noise may amplify the role of the latter against the former 
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by properly modifying the invariant measure of each map B. In fact, any 
synchronization effect is obviously absent in uniformly hyperbolic maps, e.g. 
the Bernoulli shift. 

In this paper we aim to extend this analysis to spatially extended dynam- 
ical systems under the influence of spatio-temporal noise. In particular, we 
have focused our attention on some specific coupled map lattice (CML) mod- 
els that, in the absence of noise, exhibit dynamical phases characterized by 
robust statistical properties. 

The main issues contained in this manuscript are hereafter summarized: 

• A necessary condition for obtaining synchronization in CML models is 
that the maximum Lyapunov exponent A becomes negative: this extends 
the validity of Pikovsky's criterion to spatially extended systems 0. 

• Additive spatio-temporal noise introduces effective discontinuities ip_pou- 
pled interval maps, that, for A < 0, yield "stable-chaotic" phases l!!J. In 
this case synchronization iSj-mled by a further indicator, the propagation 
velocity of information Vf 

2 Models and Tools 

The general model of CML driven by the same realization of spatio-temporal 
noise is defined by the the following two-step evolution rule: 

x\ = (1 - e)x\ + + x\ +1 ) 

xl +1 =m)+a- V l (1) 

The real state variable x\ depends on the discrete space and time indices 
i = 1,2, ■■■,L and t = 1,2, ■ ■ ■ ,T, respectively. The diffusive coupling e 
varies in the interval [0,1], a is the amplitude of the normalized space-time 
dependent random variable r]\, uniformly distributed in the interval [-1,1], and 
fix) is a mapping from S £ R onto the interval 1 GR. 

Following Ref. til, one can consider symmetric maps of the form shown in 
Fig. 1 

f(x) = tanh(Ax) - B tanh(Cx) (2) 

with A — 20, B = 1.5 and C — 2 and where S = R and T = [—a, a] 
(a ~ 0.684). The presence of the noise term of amplitude a will obviously 
enlarge the support of the invariant measure to the interval [—a — a, a + a]. 
This will be not a problem for map (||) since S is thw whole real axis. So far, 
the corresponding CML model is certainly not one of the most investigated 
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Figure 1. The map f(x) = tanh(Ar) - B tanh(Cz) for A = 20, B = 1.5 and C = 2. 



problems in the literature. In fact, CML are usually defined using maps 
/ of the interval [0,1] onto itself, like the logistic map, that are known to 
exhibit specific dynamical features, yielding many interesting statistical. and 
dynamical effects in the presence of diffusive coupling (e.g., see Ref. 
In this case X = S = [0,1] and the dynamic rule ([!]) can maintain the state 
variable x\ inside S only by adopting some further recipe: for instance x* +1 — > 



+ 1 ( xl +1 -> x* +1 - 1) if x* +1 < 0* +1 > 1). This rule modifies 



the probability measure of the state variable x\ w.r.t. the noise-free case. 
Herzel and Frcund B showed that, in low-dimensional systems, a similar recipee 
breaks the symmetry of the noise process, giving rise to an effective state- 
dependent noise, with nonvanishing mean. Later, Lai and Zhou □ have shown 
that, even maintaining symmetric noise, synchronization can be achieved, 
provided sufficiently large noise amplitudes amplify the role of contracting 
regions in phase space. 

This notwithstanding, for what concerns synchronization in high dimen- 
sional dynamical systems like CML, it is more interesting to point out that the 
above mentioned recipee introduces effective discontinuities. In fact, two state 
variables, whose values are close to one another and to one of the boundaries, 
can be driven by the addition of spatio-temporal noise at a distance of order 
1. This amounts to a strong non- linear effect, that is_highly reminiscent of 
the mechanism yielding stable chaos in CML (see Ref. LL3). Here it is induced 
by noise, while in the just quoted CML model the discontinuity was built in 
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the map /. As we are going to show in Sec. 3, the only distinctive feature 
of these scenarios is that the final attractor is a stable periodic orbit in the 
deterministic CML and a stochastic one, in the noisy model. 

In the following section we discuss synchronization in both of the above 
mentioned models of noise-driven CML. For this purpose, here it is worth 
defining the main quantitative indicators that we have used for working out 
our analysis. A first sensible quantity is the maximum Lyapunov exponent, A, 
measured according to the method outlined in Ref. u. We want to remark that 
this quantity is defined for deterministic dynamics (i.e. the noise-free case, 
(7 = 0) and it measures the average exponential expansion rate of nearby 
trajectories within a linear stability analysis. For a ^ we are faced with 
stochastic trajectories and it is not a priori obvious that A is still a well-defined 
quantity. On the other hand, the linearized evolution equations in the tangent 
space of Eq. ([[]) do not depend explicitely on r\\. Accordingly, they are formally 
equivalent for both the noisy and the noise-free cases. However, the presence 
of noise modifies the dynamical evolution of the system and therefore also the 
tangent space dynamics. Careful numerical simulations show that A remains 
a well defined asymptotic quantity, measuring the exponential expansion rate 
of infinitesimal perturbations of the stochastic evolution. In particular, its 
value is found to depend on a, but not on the realization of noise. 

Another relevant indicator for the study of spatially extended dynami- 
cal systems is the average propagation velocity of finite amplitude perturba- 
tions y 

V F = hm km ™> 

t^oo L^oo t 

where 

F(t) = maxjz : 1 < i < L; \x\ - y\\ > oj . (4) 

The trajectories x\ and y\ represent the synchronous evolution according to 
Eq. (0), with the same realization of noise 77*, of two initial conditions that 
differ at time t = for finite amounts Si ~ 0(1) only inside a space region of 
size S: y° = [x? + J,, modi] for \L/2 — i\ < S, otherwise y° = . 

In the noisy model the average in (||) has to be performed over different 
initial conditions and noise realizations. It is worth stressing that Vf is as- 
sociated to the nonlinear mechanisms of information propagation and it can 
have a finite value even for non chaotic evolution EEI. Accordingly, in order to 
have no information flow in a spatially extended system not only A but also 
Vf should vanish. 



(3) 
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Figure 2. Lyapunov exponent A as a function of the amplitude of the noise a for a chain 
of 1000 coupled maps of the type reported in Fig. [l] for e = 0.1. The A- values have been 
obtained by integrating for 130, 000 time steps. 



3 Synchronization of Replicae 

Two natural quantities can be introduced for identifying the synchronization 
of replicae: the first passage time ti(A), i.e. the time needed for two orbits, 
starting from different initial conditions xf and y®, to approach each other at 

a distance d(t) = ■h J2iLi \ x \ ~ vl\ smaller than a given threshold A (usually 
assumed much smaller thaa unity); the time 12(A) during which their distance 
remains smaller than A 113. Since, in general, these times should depend on 
the initial conditions and on the noise realizations, in what follows we shall use 
the same symbols for denoting the averaged quantities. In all the numerical 
calculations hereafter reported we have used A ~ 10~ 8 — 10 -10 . We have also 
verified that results are not affected by the choice of A, provided it is small 
enough, tipically A < 10~ 7 . 

As a first example let us consider the CML with the map / shown in Fig. 
fj In this case, we have observed synchronization for all the e-values that we 
tested. 

For instance, with a diffusive coupling value £ = 0.1 we have found that 
synchronization is obtained for a > a s = 1.44. In Fig. we show that A is a 
function of a and it becomes negative at a s . This is the first evidence of the 
validity of Pikovsky's criterion also for spatially extended systems. 

In principle, one has to consider that t\ (A) is also a function of the system 
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size L. For a < <r s , i.e. A > 0, we find that ti(A) diverges exponentially with 
the system size L : Ti(A) ~ e L ^. The scale factor £(e,a) plays the role 
of an effective correlation length and is found to be finite for any choice of 
the parameters uJ. This exponential divergence can be explained by a simple 
argument. Despite the spatial coupling, the combined effect of chaos and noise 
makes an effective number of degrees of freedom, L/£, in the replicae evolve 
independently. Accordingly, the probability that the replicae get closer than 
a distance A should be proportional to A L ^ and Ti(A) can be reasonably 
assumed to be inversely proportional to this probability. 

At variance, a logarithmic dependence of ti(A) on L is found for a > a s . 
Numerical simulations show that the number of regions made of a few synchro- 
nized sites increases as time flows. Moreover, once a synchronized region is 
formed, it grows linearly in time, until all regions merge and synchronization 
sets in over the whole lattice. Upon these observations, one can introduce 
a simple model accounting for the logarithmic dependence of t\ on L. An 
effective rate equation for the number of synchronized sites, n(t), can be con- 
structed by assigning a probability p for the formation of new synchronized 
sites and a rate 7 for the linear increase of synchronized regions: 

^= 7 + p(L-n) , (5) 

with < n(t) < L. Solving this equation with the initial condition n(0) = 0, 
one obtains an estimate of t\ by imposing the condition n(ri) = L: 



T\ = — In 
P 



pL 

— + 1 

7 



(6) 



Note that a logarithmic dependence with L (see, e.g., the inset of Fig. ^) 
is consistent with the condition ^ >> 1, that can be always achieved for 
sufficiently large values of L. 

For what Vf is concerned, in this case it is not expected to provide any 
additional information w.r.t. A. Actually, Vp is found to be positive for 
a < a s , while it vanishes for a > a s . 

Different scenarios are obtained by considering the second kind of noisy 
CML described in the previous section. Here we consider the logistic map at 
the Ulam point 

f{x) = 4z(l - x) (7) 

Let us first study the case e = 1/3: for large enough values of L, A and Vf 
are positive for any a. However, their values are found to vary with a, as well 
ri and Ti that remain finite up to a ~ 0(1). Again, t± is found to diverge 
exponentially with L. 
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Figure 3. Behaviour of the propagation velocity Vf (circles) and of the maximum Lyapunov 
exponent A (triangles) versus noise amplitude a for coupled logistic maps with e = 2/3. 
Both quantities have been computed for L = 1024, averaging over 10 3 initial conditions 
each one followed for O(10 5 ) time steps. The dashed lines are a guide for the eyes. 
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Figure 4. Exponential scaling of t\ versus L for coupled logistic maps with e = 2/3 at 
different values of the noise amplitude: a = 0.3 (empty circles), 0.32 (filled circles), 0.35 
(empty diamonds), 0.38 (filled diamonds). The inset shows the logarithmic scaling of t\ 
versus L for a = 0.45 . The values of t\ have been computed with A = 10~ 10 and averaged 
over 10 3 initial conditions. 
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Figure 5. Inverse of the scale factor l/{ as a function of cr" 1 - a' (cr c = 0.4) for CML 
logistic maps with e = 2/3 and 0.30 < a < 0.38 . 



A new interesting scenario occurs for e = 2/3, where, for L > O(10), A is 
negative for a > 0.27 (see Fig. ||). Therefore, one expects that synchronization 
occurs and this is indeed the case although, for 0.27 < a < 0.4 , Ti(A) is still 
found to increase exponentially with L (see Fig. ||). On the other hand, for 
a > 0.4 ti(A) grows logarithmically with L (see the inset of Fig. ^J. 

This transition between two completely different dynamical regimes is 
ruled by Vf, that drops to zero at a c w 0.4 (see Fig. ||) . In fact, for a < er c , 
a positive Vf implies that any perturbation of finite amplitude propagates 
through the lattice. Despite the linear mechanism is no more active (A < 0), 
we find that the non-linear mechanism of information propagation maintains 
the exponential dependence of t\ on L, thus still implying that an effective 
number of degrees of freedom evolves independently. In particular, we find 
numerical evidence of the relation oc (<r _1 — a^ 1 ) for a < a c (see Fig. ||); 
this indicates that the scale factor £ diverges at the transition point. 

Eventually, only when d(t) becomes sufficiently small the linear mecha- 
nism acts as a stabilizing factor on the dynamics of the replicae and synchro- 
nization is achieved. 

For a > 0.4, where also the nonlinear mechanism is absent, numerical 
simulations exhibit the same scenario described in the "symmetric map" CML 
for a > a s . 

One can counclude that in this last case the conditions A < and Vf — 
have to be fulfilled independently for guaranteeing noise-driven synchroniza- 
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tion of two replicae within a "reasonable" time span in a CML of large but 
finite size L. 

We want to point out that this scenario is not peculiar of the logistic map 
for e — 2/3. Actually, we have verified that it holds in a finite interval of 
values around e = 2/3. Moreover, a very similar situation is pbtained when 



where b — 2.7, a — 0.07 and c = 0.1. This case is particularly interesting for 
a twofold reason. It has been shown that the noise-free case has a negative 
A for any value of e and, moreover, it exhibits a peculiar transition at e c ~ 
0.6 from_a frozen disordered phase with Vf = to a chaotic phase with 
Vf > Ej. At variance with the case of coupled logistic maps, when noise is 
added A remains negative for any value of a. The phase transition disappears, 
because, independently of e, Vf is found to be positive for very small noise 
amplitudes. On the other hand, by increasing a up to a critical value u c (e) , 
Vf is found to drop again to zero not only below, but also above e c : for 
instance, one has a c « 0.16 for e = 0.58 and a c « 0.18 for e — 0.62. In both 
of these cases we have recovered the same kind of mechanisms characterizing 
the synchronization transition discussed for the case of coupled logistic maps 
with e = 2/3 (see Fig. ^). In this sense, it seems reasonable to conjecture that 
this phenomenon is present in a wide class of spatially extended dynamical 
systems and that it can be characterized making use of the same kind of 
analysis discussed in this paper for CML models. 

As a final remark, we would like to stress that all this scenario is rem- 
iniscent of the phenomenology associated to stable chaos in CML, that was 
identified for model (|h li!J. Also in that case unpredictability was character- 
ized by transients increasing exponentially with L when Vp > 0. The analogy 
is strenghtened by the fact that additive noise induces "jumps" in the de- 
terministic dynamics of the system, similar to the discontinuities present in 
model (JsJ) . As a matter of fact, the presence of strong variations in the slope 
of the local p&p is fundamental in order to observe non-linear propagation 
mechanisms 113. 
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